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ABSTRACT 

Gravitational lensing by massive galaxy clusters is a powerful tool for the discovery and study of high red- 
shift galaxies, including those at z>6 likely responsible for cosmic re-ionization. Pello et al. recently used 
this technique to discover a candidate gravitationally magnified galaxy at z=10 behind the massive cluster lens 
Abell 1835 (z=0.25). We present new Keck (LRIS) and Spitzer Space Telescope (IRAC) observations of the 
z=10 candidate (hereafter#1916) together with a re-analysis of archival optical and near-infrared imaging from 
the Hubble Space Telescope and VLT respectively. Our analysis therefore extends from the atmospheric cut-off 
at A bs— 0.35/im out to A b s — 5/zm with Spitzer/IRAC. The z=10 galaxy is not detected in any of these data, 
including an independent reduction of Pello et al.'s discovery H- and /T-band imaging. We conclude that there 
is no statistically reliable evidence for the existence of #1916. We also assess the implications of our results 
for ground-based near-infrared searches for gravitationally magnified galaxies at z>7. The broad conclusion is 
that such experiments remain feasible, assuming that space-based optical and mid-infrared imaging are avail- 
able to break the degeneracy with low redshift interlopers (e.g. z^2-3) when fitting spectral templates to the 
photometric data. 

Subject headings: cosmology:observations — early universe — galaxies: evolution — galaxies: formation — 
infrared: galaxies 



1. INTRODUCTION 

Observations of distant QSOs (Becker et al. 2001; Fan et 
al. 2002) and the cosmic microwave background (Kogut et al. 
2003) together suggest that the universe was re-ionized some- 
where between z— 6 and z— 20. Searching for the sources of 
re-ionizing photons is currently an intense observational ef- 
fort. Most searches naturally concentrate on luminous sys- 
tems, i.e. QSOs and luminous galaxies, at z~6-8 as these 
should be easier to detect than less luminous and more distant 
objects. However QSOs likely produced insufficient photons 
to accomplish re-ionization alone (Fan et al. 2001; Barger et 
al. 2003), and the same may be true of luminous (L;>0.3L* =3 8 ) 
galaxies based on small samples from the Hubble Space Tele- 
scope Ultra Deep Field (hereafter HST UDF; Bouwens et al. 
2004; Bunker et al. 2004; Yan et al. 2004). This raises the im- 
portant possibilities that re-ionization either occurred much 
earlier, or that the bulk of the re-ionizing photons were emit- 
ted by sub-luminous galaxies, i.e. L<().IL*. 

The UDF studies operate close to the detection threshold 
of the deepest optical/near-infrared imaging available. It is 
therefore difficult to envisage substantial progress in the de- 
tection of more remote and/or less luminous galaxies via deep 
imaging of "blank fields" with the current generation of tele- 
scopes. With the advent of the James Webb Space Telescope 
still some years ahead, the magnifying power of massive 
galaxy cluster lenses is therefore a much needed boost for the 
discovery power of HST and large ground-based telescopes. 
Indeed, the galaxy redshift record has been broken on several 
occasions with the help of the gravitational magnification of 
distant galaxies by foreground galaxy clusters (Mellier et al. 
1991; Franx et al. 1997; Hu et al. 2002; Kneib et al. 2004). 
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The faint end of the luminosity function of Lyman-a emitters 
at z=5 has also been constrained with the help of gravitational 
lensing (Santos et al. 2004; Ellis et al. 2001). Extension of 
these techniques to z;>7 is therefore an important element of 
observational studies of cosmic re-ionization . 

Pello et al. (2004 - hereafter P04) reported a gravitation- 
ally magnified (/i^25-100) galaxy at z=10 (hereafter #1916, 
following P04's nomenclature) behind the foreground galaxy 
cluster A 1835 (z=0.25). This interpretation is based on non- 
detection in optical imaging from the ground (3cr limits in 
a 0.6" diameter aperture: V>27.4, R>21.5, 7>26.9) and 
space (3(7 limit in a 0.2" diameter aperture: 7?702>27.2), 
and the shape of the continuum at A bs>l/im (using a 1.5" 
aperture: (J-H)>0.6, (H-K)=-0.5±0A) which is reminis- 
cent of the Lyman-break selection technique (Steidel et al. 
1996). P04 corroborated the putative Lyman-break redshifted 
to A b s ~1.3/im with an emission line at A b s =1.3375/im with 
integrated flux of (4.1±0.5)xl0" 18 ergcm" 2 s _1 , which they 
interpret as Lyman-a. Lower redshift interpretations of the 
line ([Oil] at z=2.59; [OIII] at z=1.68; Ha at z=1.04) were 
discarded by P04 largely on the basis of the low probability of 
solutions at z<^7 when fitting synthetic spectral energy distri- 
butions (SEDs) to their photometric data. The most likely of 
these lower-redshift solutions (z=2.59) was further excluded 
on the basis of the dust extinction required to fit the photo- 
metric data (Ay>2), and the absence of doublet structure in 
the observed emission line. 

Based solely on the photometry (optical non-detection, red 
(J-H) and blue (H-K) colors) the z=10 interpretation of 
#1916 is plausible. However P04's preference for this so- 
lution over the lower redshift alternatives was controversial 
from the outset. For example, the emission line does not have 
the characteristic P-Cygni profile of Lyman-a, and the dif- 
ferent photometric apertures adopted in the optical (0.6" - 
smaller than the ground-based seeing disc) and near-infrared 
(1.5" - 3x the seeing disk) may suppress the likelihood of 
lower-redshift solutions when fitting synthetic SEDs. Bremer 
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et al. have also suggested that #1916 may either not exist, or 
be intrinsically variable, based on their non-detection in the 
//-band with NIRI on Gemini-North, H(3a)>26, in contrast 
to P04's 4cr detection of i/=25.00±0.25 with ISAAC on VLT. 
The spectroscopic identification of #1916 is also in doubt. 
Weatherley et al. (2004) re-analyzed P04's spectroscopic data 
and failed to detect the emission line at A bs=L3375/im, cit- 
ing spurious positive flux arising from variable hot pixels in 
the ISAAC array as the likely source of the discrepancy. 

A 1835 has been used previously as a gravitational tele- 
scope, for example targeting sub-millimeter galaxies (here- 
after SMGs; e.g. Smail et al. 2002) and galaxies with 
extremely red optical/near-infrared colors (Smith et al. 
2002). One of the galaxies detected in these surveys, 
SMMJ 14011+0252, lies at z=2.56 and suffers an estimated 
extinction of 1.8<^4y<6.5 (Ivison et al. 2000). Bearing in 
mind recent discovery of galaxy groups at z^2-3 associated 
with gravitationally lensed SMGs (Kneib et al. 2004; Bo- 
rys et al. 2004), and the strong clustering of SMGs (Blain 
et al. 2004), #1916 is plausibly at a similar redshift to 
SMMJ 1401 1+0252, and may also be obscured by dust. Fur- 
ther circumstantial evidence for a lower redshift interpretation 
of #1916 comes from Richard et al. (2003) who used the same 
spectroscopic data as presented by P04 to discover a strongly 
reddened star-forming galaxy at z=1.68. This redshift coin- 
cides with the [OIII] interpretation of P04's putative emission 
line at A b s =L3375/im. 

In this paper, we address three questions: (i) is #1916 
at z^2-3?; (ii) is #1916 intrinsically variable?; (iii) does 
#1916 exist? These tests exploit new optical and mid- 
infrared observations using the Keck-I 10-m telescope and 
the Spitzer Space Telescope respectively, plus an independent 
reduction of P04's H- and /^-band imaging data from VLT. 
Throughout this article we assume that the emission line at 
A bs=1.3375yum is a false detection (Weatherley et al. 2004). 
In ^2] we present the new Keck and Spitzer data, explain in 
detail the re-reduction of the archival VLT/ISAAC data, and 
summarize the archival Hubble Space Telescope (HST) data. 
Then in we describe the analysis and key results, focus- 
ing on the three questions posed above; this section closes 
with a summary of the current observational status of #1916. 
Finally, we discuss the implications of our results for future 
ground-based near-infrared searches for galaxies at z>7 (S|4} 
and summarize our conclusions in (|5J. 

We assume //o=65kms" 1 Mpc" 1 , £!m=0.3 and fl\=0.1 
throughout. Unless otherwise stated all error bars are at 
la significance; photometric detection limits are at 3a sig- 
nificance; upper and lower limits on colors are based on 
3a detection thresholds in the non-detection filter. Mag- 
nitudes are stated in the AB system; conversion between 
the AB and Vega systems for the specific filters used in 
this paper are as follows:- A£=BAB-Bveea= - 0.1, Ay=0.1, 
A R =0.2, A F702W =0.3, A 7 =0.5, A y =0.9, A H =1A, A K =1.9, 

A3.6 / im=2.8, A4.5 /m =3.2. 

2. OBSERVATIONS AND DATA REDUCTION 

We describe new and archival observations of #1916 in 
order of increasing wavelength, spanning the observed opti- 
cal, near-infrared and mid-infrared: spectroscopy with LRIS 
on Keck-I (0.35<A obs <0.95/im); imaging with WFPC2 on- 
board HST (A o bs=0.7^m); near-infrared imaging with ISAAC 
on ESO's VLT (A ob s=L6^m and 2.2^m); IRAC/Spitzer ob- 
servations at A bs=3.6/im and 4.5/im. 

The detection of any optical flux from #1916 would elim- 



TABLE 1 
Summary of Photometry 



Filter Telescope/ 
Instrument 



FWHM 

(") 



Aperture Photometry" 
Pello et al.* This Paper 



V CFHT/12k 0.76 

R CFHT/12k 0.69 

R-,02 #S77WFPC2 0.17 

/ CFHT/12k 0.78 

J VLT/ISAAC 0.65 

H VLT/ISAAC 0.50 

K VLT/ISAAC 0.38 

3.6/im SpitzermAC 1.7 

4.5//m Spitzer/IRAC 1.7 



> 27.5 (0.6") 
>27.6(0.6") 

>27.2(0.2") c >27.0(0.5") 
> 26.0 (0.6") 
>25.6(1.5") 
25.00±0.25(1.5") > 25.0(1.5") 
25.51±0.36(1.5") > 25.0(1.5") 
>24.3(5.1") 
>24.3(5.1") 



a Each number in parentheses is the diameter of the aperture used for the 
respective photometric measurements. In il4lP04's optical non-detections are 
re-scaled to a photometric aperture of 2" diameter (~3 times the seeing disc): 
V>26.2, S>26.3, />24.7. 

b We convert all of P04's optical detection limits and their lcr 7-band detec- 
tion to 3cr limits. 

c P04 do not state whether their R702 detection limit is in the Vega or AB 
system. We have assumed the former and converted it to AB in this table. 
P04 also do not explain how they reduced the WFPC2 data, specifically 
whether the final pixel scale was different from the native 0.0996" /pix of 
the WFC detectors. In this table we have assumed that the four pixels over 
which this detection limit is measured (see P04 §2.1) subtend a solid an- 
gle of 0.0996" X 0.0996". Given these uncertainties and the absence of this 
detection limit from P04's Table 1 and Fig. 3, we ignore P04's limit when 
attempting to reproduce their photometric redshift results in ij4] 



inate the z=10 interpretation (e.g. Stern et al. 2000). In con- 
trast, optical non-detection would have several alternative in- 
terpretations, including: a galaxy at z=10 as per P04; a dusty 
galaxy at z^2- 3; n on-existence of #1916. The mid-infrared 
observations ( 32.4> should therefore help to constrain the 
amount of energy re-radiated by dust in the z~2-3 interpre- 
tation, and our re-reduction of P04's VLT data also helps to 
clarify the possibility that #1916 may not exist, or may be 
variable (Bremer et al. 2004). The Keck and Spitzer data de- 
scribed in this section were collected and analyzed in parallel 
with those presented by Bremer et al. (2004). 

2.1. Keck Spectroscopy 

As part of a broad effort to secure spectroscopic redshifts 
of gravitational arcs spanning several observational programs 
(Smith et al. 2001, 2002, 2005; Sand et al. 2002, 2004, 2005; 
Edge et al. 2003; Sharon et al. 2005), we observed A 1835 
with the Low Resolution Imager Spectrograph (LRIS; Oke et 
al. 1995) in multi-slit mode on the Keck-I 10-m telescope 4 on 
UT 2004 March 29. A single mask was observed, including a 
slit targeting #1916. The purpose of this slit was to search for 
line emission in the 0.35<A o b s <0.95/im wavelength range. 
For example, if #1916 does indeed lie at z~2.6 (JD, then 
Lyman-a may be detectable at A o b s ~0.44/im. 

The observations totaled 3.6-ks, split into two exposures, 
using the D560 dichroic with the 400/8500 grating and the 
400/3400 grism. On the red side the spectral dispersion was 
1.86A/pixel with a pixel scale of 0.214"/pixel and on the 
blue side the spectral dispersion was 1.09 A/pixel with a pixel 
scale of 0.135"/pixel. Overhead conditions were moderate 
(FWHM^l"), and probably not photometric, however a flux 

4 The W. M. Keck Observatory is operated as a scientific partnership 
among the California Institute of Technology, the University of California, 
and NASA. 
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FIG. 1 . — Infrared imag es of the z=10 candidate at 1 .6, 2.2, and 3.6/im respectively. The two left panels are based on our independent re-reduction of P04's 
ISAAC data described in j|2.3l The white circles mark the position of #1916 from P04 - there is no obvious sign of flux in any of these panels. Formal 3a 
detection limits are: H>25, K> 25, F<0J5p]y. North is up and East left. Each panel is 23" X 16". 



calibration was obtained using the spectrophotometric stan- 
dard star HZ44 (Oke et al. 1990). The data were de-biased, 
flat-fielded, sky-subtracted, extracted and calibrated in a stan- 
dard manner within IRAF 5 . 

No flux at observed optical wavelengths has yet been de- 
tected at the position of #1916 (P04; Lehnert et al. 2004; 
32.2I >. We therefore did not expect to detect any contin- 
uum emission, and concentrated instead on searching for faint 
emission lines of large equivalent width. Visual inspection of 
the reduced 2D data revealed neither continuum nor line emis- 
sion. To estimate the sensitivity limit we extracted a ID trace 
from the center of the slit corresponding to the full width of 
the seeing disk, and estimated the 3a detection limit per 5-A 
spectral resolution element to be ~4.5 x 10~ 19 ergs _1 cm" 2 at 
A o b s =0.44/im - i.e. the observed wavelength at which Lyman- 
a would be found if #1916 is at z=2.6. 

2,2. Archival Hubble Space Telescope Imaging 

A 1835 has also been observed through the F702W filter 
with the WFPC2 camera on-board HST b . We refer the reader 
to Smith et al. (2005) for details of these data and their reduc- 
tion. P04 do not detect #1916 in these data (Table[T}, although 
they neither explain how they reduced the data nor how the de- 
tection limit was calculated. Here, we use Smith et al.'s (2005) 
reduced frame which has a pixel-scale of 0.0498"/pixel after 
drizzling. To simplify the analysis, we re-bin the data back 
to the original pixel-scale of 0.0996"/pixel to minimize the 
impact of pixel-to-pixel correlations in the background noise 
when estimating the sensitivity limit of the data. 

Visual inspection reveals no obvious flux at the position 
of #1916 in these data. To quantify this non-detection, we 
follow the same procedure as P04 - we measured the back- 
ground noise in apertures placed randomly into blank sky re- 
gions of both frames near to the position of #1916. We ensure 
consistency with other wavelengths by matching the diameter 
of these apertures to a diameter 3 x that of the seeing disk, 
i.e. 0.5". We obtain a 3a sensitivity limit in that aperture of 
#702=27.0 (TableHJ. 

2.3. Archival VLT/ISAAC Imaging 

5 IRAF is distributed by the National Optical Observatories, which is op- 
erated by the Association of Universities for Research in Astronomy, Inc. 
(AURA) under cooperative agreement with the National Science Foundation. 

6 Based in part on observations with the NASA/ESA Hubble Space Tele- 
scope obtained at the Space Telescope Science Institute, which is operated 
by the Association of Universities for Research in Astronomy, Imc, under 
NASA contract NAS 5-26555. 



Deep near-infrared imaging of A 1835 was obtained by 
P04 with the ISAAC 1024x1024 Hawaii Rockwell array on 
ESO's 8-m VLT 7 in February 2003. We reduced indepen- 
dently the H- and /^-band data using standard IRAF tasks, 
paying careful attention to rejection of cosmetic defects in 
the ISAAC array, including bad pixels, and to conserving the 
noise properties of the data. We detail key features of the data 
reduction below. 

(i) Flat-fielding and sky-subtraction were combined into a 
single step using a median of the eight frames temporally 
adjacent to each science frame. We refer to this step as 
"flat-fielding", and to the rolling temporal median frames 
as "sky-flats". 

(ii) Flat-fielding was performed twice. First on the dark- 
subtracted frames which were then registered and aver- 
aged to produce a first-pass reduced frame. This frame 
was then used to mask out flux from identified sources 
from the individual dark-subtracted frames, and these 
masked frames were then used to construct the sky-flats 
in a second-pass reduction. This approach minimizes 
the loss of flux from objects with angular extents com- 
parable with the size of the dither pattern. This is im- 
portant when searching for faint objects along lines-of- 
sight through the crowded cores of rich galaxy clusters 
because the light from bright cluster galaxies effectively 
form a spatially varying background against which the 
faint sources are detected. The goal of the second-pass 
flat-field is to conserve this "background". 

(iii) Independent bad pixel masks were made by sigma- 
clipping both the darks and the sky-flats. The former 
identifies 22, 020 pixels (2. 1 % of the total array) as static, 
i.e. "bad", and the latter identifies the same 22,020 pix- 
els plus an additional 12,373 pixels (a further 1.2% of 
the total array) as bad. The latter mask was adopted as 
the fiducial bad pixel mask. 

(iv) Detector bias residuals were removed by subtracting the 
median along rows in individual flat-fielded frames after 
masking identified sources, in a manner similar to that 
described by Labbe et al. (2003). 

(v) The individual frames were integer pixel aligned. This 
has the important benefit of minimizing pixel-to-pixel 

7 Based in part on observations collected with the ESO VLT-UT1 Antu 
Telescope. 
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correlations in the noise properties of each frame and 
thus the final stacked frame. Calculation of the back- 
ground noise is therefore simplified relative to a reduc- 
tion scheme based on sub-integer pixel alignment of the 
individual frames. 

(vi) No frames were rejected when making the final combi- 
nation of the reduced, aligned frames. Two versions of 
the final stacked frame were made: a straight average and 
a weighted average - the weight of an individual frame 
was proportional to (a xrms)" 2 , where a is the FWHM 
of the seeing disk, and rms is the root mean square per 
pixel of the noise in each frame. The weighted version of 
the final frame has slightly better image quality than the 
straight average. We therefore adopt it for the analysis 
described below. 

The final reduced H- and /f-band frames have seeing 
of FWHM=(0.45±0.01)" and FWHM=(0.34±0.02)" respec- 
tively. Photometric calibration was achieved with the standard 
star observations that were interspersed with the science ob- 
servations as part of P04's original program. We show extracts 
from the reduced H- or /f-band frames in Fig. ^ Visual in- 
spection of both the fits frames and Fig.H revea l s no obvious 
flux at the location of #1916 in either the H- or /f-band. Fol- 
lowing P04 we again randomly insert 1 .5" diameter apertures 
(roughly 3 times the seeing disc) into blank sky regions near 
to the position of #1916, obtaining 3a sensitivity limits in this 
aperture of: H=25 .0 and £=25 .0. 

2.4. Spitzer///MC Imaging 

A 1835 was observed with the InfraRed Array Camera 
(IRAC; Fazio et al. 2004) on-board Spitzer* on UT 2004 Jan- 
uary 16 in the 3.6, 4.5, 5.8 and 8.0^m channels. Here we 
discuss the two shortest wavelength, more sensitive, observa- 
tions. Twelve and eighteen 200-second exposures were ac- 
cumulated at 3.6/im and 4.5/im respectively, using the small- 
step cycling dither pattern. The Basic Calibrated Data (BCD) 
were combined using custom routines to produce the final 
stacked frame with a pixel scale of 0.6"/pixel. We re-binned 
the data back to the original pixel scale of 1 .2" /pixel to elim- 
inate correlations in the background noise. 

Visual inspection of the final frames again indicates that 
there is no flux at the position of #1916 (Fig.Q. To quantify 
this non-dete ction we follow the same procedure as P04, as 
described in SI2.2I We used 5.1" diameter apertures, i.e. 3x 
the seeing disk of the IRAC observations to obtain 3a sensi- 
tivity limits of: F(3.6/im)=0.75/iJy and F(4.5//m)=0.75/iJy 
respectively. 

3. ANALYSIS AND RESULTS 

The objective of this section is to answer the three questions 
posed in ffl (i) is #1916 at z~2-3?; (ii) is #1916 intrinsically 
variable?; (iii) does #1916 exist? Preliminary inspection of 
the data in ^indicates that no flux is detected at the posi- 
tion of #1916 at any wavelength to date. Combining this with 
Bremer et al.'s more sensitive non-detection of //(3ct)>26.0, 
it is tempting to leap to the third question and reply "no". We 
adopt a more conservative approach. 

8 This work is based in part on observations made with the Spitzer Space 
Telescope, which is operated by the Jet Propulsion Laboratory, California 
Institute of Technology under a contract with NASA. 



3.1. Is #1916 at z~2-3? 

This test concentrates on the optical data because the de- 
tection of any flux shortward of the putative Lyman limit of 
a galaxy at z— 10 would immediately discount that interpre- 
tation. The red observed optical/near-infrared spectral en- 
ergy distribution described by P04 could then be naturally ex- 
plained by a dusty galaxy at z~2-3, perhaps associated with 
the SMGs that lie within -30" (~200-300kpc in projection 
at z~2-3) of #1916 (Ivison et al. 2000; Smail et al. 2005). 

Our new non-detection of #1916 with LRIS ( SI2.1> . coupled 
with confirmation of P04's non-detection with HST/WFPC2 
and Lehnert et al.'s non-detection in the V-band with 
VLT/FORS are mutually consistent in the sense that no optical 
flux has been detected at this position to date. However these 
non-detections are consistent with all of the following: z=10, 
extreme dust obscuration at z^2-3, an intrinsically variable 
source, and non-existence. The result of this test is therefore 
inconclusive. 

3.2. Is #1916 Intrinsically Variable? 

The objective of this section is to test Bremer et al.'s (2004) 
proposal that #1916 is intrinsically variable. If P04's photom- 
etry (//=25.00±0.25 and K=25. 51 ±0.51) is reproducible us- 
ing our independent reduction of their near-infrared data, then 
the variable hypothesis would be supported. If not, t hen t he 
idea that #1916 does not exist would gain credibility ( 33.3> . 

We attempt to reproduce P04's analysis using SExtractor 
(Bertin & Arnouts 1996). SExtractor was configured to lo- 
cate all sources with at least 7 pixels that are >0.75<r per 
pixel above the background - i.e. a signal-to-noise ratio of 
>2 per resolution element, based on the //-band seeing disk 
of FWHM=0.45±0.01" ( $T& and the 0.15"/pix scale of the 
ISAAC pixels. We also smoothed the data with a gaussian fil- 
ter that matched the FWHM of the observed point sources, i.e. 
a gaussian of FWHM=3 pixels. In this configuration SExtrac- 
tor failed to detect a source at the position of #1916. We there- 
fore experimented with different smoothing schemes, both in- 
creasing and decreasing the full width of the gaussian filter. A 
"detection" was only possible with the smallest available fil- 
ter - FWHM=1.5 pixels, i.e. half the width of the seeing disk 
- yielding //=25.3±0.6. Experimentation with block filters 
produced similar results in that a "detection" was not possible 
with any of the standard SExtractor block filters: 3x3, 5x5, 
7x7 pixels. We also analyzed the ZT-band data in exactly the 
same manner and failed to detect anything at the position of 
#1916 with any gaussian or block filter. 

The //-band segmentation map produced when smoothing 
with the FWHM=1.5 pixel gaussian reveals that the "detec- 
tion" is very elongated, with a width of 1-2 pixels and a length 
of ^5 pixels. The orientation of these pixels is consistent with 
the orientation of #1916 reported by P04. It is important to 
stress that the motivation for filtering data with a kernel that 
matches the resolution element of the data is to suppress false 
detections. The collection of pixels identified by SExtractor 
at the position of #1916 was only "detectable" with a smooth- 
ing kernel that has a linear scale half that of the resolution 
element of the data. It is therefore instructive to consider how 
many such ^2a blobs exist within the ISAAC data. In a single 
1.5" diameter aperture (i.e. matching that used for the pho- 
tometry described above) placed randomly in these //-band 
data, there is a 5% chance of detecting a 2a noise fluctua- 
tion - i.e. a spurious detection. However the ISAAC array 
(1024xl024pixels, each pixel 0.15"x0.15") contains of or- 
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der 10 independent photometric apertures of 1.5" diameter. 
The //-band frame therefore contains ~500 noise fluctuations 
of 2c significance. 

Sadly, the only reasonable conclusion to draw from this 
analysis is that #1916 is not detected in our independent re- 
duction of P04's data. We therefore place 3* 7 lim its on the 
flux at this position of: //>25.0 and K>25.0 (JSJ. The only 
wavelength at which two directly comparable observations are 
available is in the //-band. Combining our non-detection with 
that of Bremer et al. (2004), we conclude that there is no ev- 
idence for variability of #1916, and that (if it exists) its H- 
band flux is fainter than H=26 at 3<r significance (Bremer et 
al. 2004). 

3.3. Does #1916 exist? 

The results of the preceding two sections were derived from 
non-detection of #1916 across the broadest wavelength range 
to date: 0.35<A o b s <5/im. We now combine all of these non- 
detections to address the question of whether #1916 exists. 
The data force us to conclude that there is no statistically 
sound evidence that #1916 exists. The balance of probability 
is that #1916 was a false detection in P04's discovery obser- 
vations. New observational data yielding statistically sound 
detections are required before P04's claim that #1916 is the 
most distant galaxy yet discovered may be resurrected. We 
consider this unlikely, but we hope to be surprised by Pello et 
al.'s forthcoming HST observations. 

4. DISCUSSION 

4. 1 . The Near-infrared Non-detection 

We first discuss possible reasons for the difference between 
our non-detection of #1916 and P04's ~3-4er near- infra red 
detections. We single out two data reduction steps ( 32.3> for 
discussion - the efficiency of bad pixel rejection and the con- 
servation of noise properties. 

4.1.1. Efficiency of bad pixel rejection 

The efficiency with which bad pixels are identified can af- 
fect estimates of the signal coming from faint sources - if 
some bad pixels are not identifie d, th ey could enhance the 
flux detected by SExtractor. In 32.31 we used two different 
methods to identify bad pixels, finding a small but potentially 
important difference between the two methods. To assess the 
impact of reduced efficiency of bad pixel identification we 
made a mask image that contained the 12,373 bad pixels iden- 
tified only in the bad pixel mask generated from the sky-flats. 
We then made one copy of the mask per science frame and 
integer-pixel shifted them to match the observed dither pat- 
tern. Finally we took the weighted average of these aligned 
mask frames and summed the pixel counts in a 1.5" diameter 
aperture centered on #1916. From this we concluded that 7 
bad pixels that are only identified in our "sky-flat" bad pixel 
maps fall within the final photometric aperture. We estimate 
that if not identified and excluded from the analysis, these 
pixels could increase the flux estimates by several tenths of a 
magnitude. 

4.1.2. Conservation of noise properties 

The approach to re-sampling (or not) the individual frames 
during the data reduction process, especially the alignment 
of individual frames affects the noise properties of the fi- 
nal stacked frame. If the individual frames are re-sampled, 
for example by sub-integer pixel aligning them immediately 



prior to producing the final stacked frame, then the noise in 
the stacked frame is correlated. Such pixel-to-pixel correla- 
tions are generally absent from integer-pixel aligned data, thus 
simplifying the noise properties of the final frame. Neither 
sub-integer nor integer pixel alignment is intrinsically cor- 
rect. The relevant issue is correct measurement of the noise 
in each case - this is critical to assess accurately the statistical 
significance of sources detected close to the sensitivity limit 
of the data. Specifically, if the pixel-to-pixel correlations in 
sub-integer pixel aligned data are not included in the error 
analysis, then the noise is under-estimated and the statistical 
significance of the detection over-estimated (Casertano e t al. 
2000). We integer pixel aligned the individual frames in 32.31 
in order to simplify the error analysis. We now estimate by 
how much we would have under-estimated the noise if we 
had sub-integer pixel aligned the individual frames and then 
ignored the pixel-to-pixel correlations when calculating noise 
level. This is achieved by simply sub-integer aligning the in- 
dividual frames and re-combining them using a weighted av- 
erage. Ignoring any resulting pixel-to-pixel correlations, we 
obtain a 3<7 threshold of Z/=25.2, which is slightly fainter than 
our threshold of Z/=25.0. 

In summary, it is plausible that the difference between our 
near-infrared non-detection and P04's detection of #1916 us- 
ing the same raw data can at least in part be explained by the 
efficiency of bad pixel identification and treatment of corre- 
lated noise. 

4.2. Implications for Future Work 

We now discuss the implications of our results for future 
work, focusing on the feasibility of searches for gravitation- 
ally magnified stellar systems at z;>7 using ground-based 
near-infrared data. We begin by noting that, had their pho- 
tometry been reliable, P04's original z=10 interpretation of 
#1916 was plausible, based solely on the photometric data. 
We therefore adopt P04's optical and near-infrared photom- 
etry as being representative of what may be expected from 
similar future experiments - i.e. optical non-detections in sev- 
eral filters based on a few hours of observations with a 4-m 
class telescope, a non-detection in a single red optical filter 
with HST, ^3-4(7 detections in two near-infrared filters, and 
possibly observations with SpitzerflRAC. Specifically, we in- 
vestigate the degeneracy between z^>7 interpretations of such 
datasets with lower redshift alternatives, and how such degen- 
eracies might be broken. 

We use Version 1.1 of HYPER-Z 9 (Bolzonella et al. 2000) 
to fit standard Bruzual & Chariot (1993) single stellar pop- 
ulation models (Burst, E, Sa, Sc, Im) to the photometric 
data. We assume a Calzetti et al. (2000) extinction law, al- 
low dust extinction in #1916 to lie in the range 0<Ay<4, 
adopt E(B-V) = 0.03 for extinction within the Milky Way 
(Schlegel et al. 1998), and use Madau's (1995) prescription 
for absorption by the inter-galactic medium. In each case 
described below, if we fit the models to all available photo- 
metric information across the full redshift range (0<z<ll), 
then we obtain x 2 <^l (all \ 2 values are quoted per degree of 
freedom) for all redshifts beyond z—1- This is because the 
models are defined to have zero flux short-ward of the Lyman 
limit (A rest =0.0912/im). At z;>6, the models also have neg- 
ligible flux short-ward of Lyman-a (A re . s t=0.1215/im) due to 
the Lyman forests and Gunn Peterson absorption. These two 

9 Available from http://webast.ast.obs-mip.fr/hyperz 
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FIG. 2. — LEFT: Reduced \ 2 as a function of photometric redshift and the filters used in the model fits. The red dashed curve has been offset +0.2 in the 
F-direction for clarity. This panel demonstrates that P04's discovery photometry (VRIJHK) is degenerate between z—2.6 and z=10 (blue solid curve). Adding 
in sensitive space-based detection thresholds from HST and Spitzer/IRAC (red dashed and black dot-dashed curves respectively) lifts this degeneracy. RIGHT: 
The best-fit spectral templates to the R-jqiJ UK-band photometry at z=2.6 (red solid) and z=10 (blue dashed). The IRAC detection thresholds are also marked to 
illustrate the power of these data to discriminate between low- and high-redshift interpretations of the shorter wavelength data. 



spectral features are progressively redshifted long-ward of the 
observed optical filters at high redshift, rendering the shorter 
wavelength detection limits irrelevant to the fits. The good- 
ness of fit is therefore systematically over-estimated at high 
redshift if all photometric information is included in the fit 
across the full redshift range. We therefore fit the models to 
the data in a series of redshift chunks - only the photometric 
information that lies longward of red-shifted Lyman-a is con- 
sidered in each redshift chunk. The results described below 
are insensitive to whether the Lyman limit is substituted for 
Lyman-a. 

First, we fit the models to P04's VRIJHK-band ground- 
based photometry (i.e. the V RI J-band detection limits and the 
HK-band detections listed in column 4 of TableQ. Note that 
P04 used different sized apertures for the near-infrared (1 .5") 
and optical (0.6") photometry respectively. We find two ac- 
ceptable solutions: z—2.6 and z— 10 (Fig. [5), neither of which 
require any dust-obscuration within #1916, the latter having 
the lower formal x 2 value. We also scale P04's ground- 
based optical non-detections to that appropriate for a consis- 
tent photometric aperture of three times the seeing disk at all 
wavelengths (V >26.2, R>26.3, />24.7). This marginally im- 
proves the goodness of fit of the z— 2.6 solution, and is oth- 
erwise indistinguishable from the original fit. We retain the 
matched photometric apertures for the remainder of this anal- 
ysis. 

To improve the discrimination between low- and high- 
redshift solutions, we add the sensitive //S77WFPC2 de- 
tection threshold of /?702>27.0 (Tabled to the photomet- 
ric constraints. We fit the synthetic SEDs to the combined 
VRRjoiUHK dataset. The goodness of fit of the z—2.59 so- 
lution is marginally worse relative to the VRIJHK-band anal- 
ysis, because the most stringent optical non-detection comes 
from the /?702-band. However, in general, the \ 2 as a func- 
tion of redshift is indistinguishable between this fit and the 
VRIJHK-band fits. This is probably because the \ 2 is dom- 
inated by the non-detections in six out of the eight observed 
filters. To test this we re-fit the models, limiting the data to 
just the R-ioiJHK-bands, i.e. the most sensitive non-detection 
(Rjqi), the reddest non-detection (J) and the two detections 
(HK). The impact of the sensitive detection limit from the 



HST data is now clearly evident in the new \ 2 distribution, as 
shown in Fig. [21 The only acceptable fit to the RjoiJHK-band 
data is at z~10. 

We show the best fit SEDs at z=2.6 and z=10, based on 
the Rj^JHK-band data in the right panel of Fig. |2] This 
demonstrates the potential power of IRAC photometry to fur- 
ther discriminate between low and high-redshift interpreta- 
tions of candidate high-redshift galaxies. The SED of low 
redshift solutions (e.g. z— 2.6) is red at A b s ^3-5/j,m, and the 
SED of high redshift solutions (z~ 1 0) is blue at similar wave- 
lengths. We add the detection limits obtained at A b s =3.6/im 
and 4.5/im to the R^JHK-band data and re-fit the spectral 
templates. The result is shown in the left panel of Fig.|5]- the 
reduced \~ of the z~2.6 solution is now in excess of 2, and 
the only redshifts which achieve an acceptable fit to the data 
are 9.5£z^l 1.5. 

In summary, a single sensitive non-detection derived 
from HST imaging, combined with deep ground-based 
near-infrared imaging and ^2,400 second integrations with 
Spitzer/IRAC at 3.6 and 4.5/im can break the degeneracy be- 
tween low and high-redshift interpretations of candidate z;>7 
stellar systems (see also Egami et al. 2004; Eyles et al. 2005). 
Whilst this is not an exhaustive study, it demonstrates that de- 
spite the demise of #1916, searches for gravitationally magni- 
fied galaxies at extremely high redshifts using ground-based 
near-infrared data remain feasible when combined with sen- 
sitive space-based optical and mid-infrared data. 

5. CONCLUSIONS 

We have analyzed new and archival observations of the 
z=10 galaxy #1916 behind the foreground galaxy clus- 
ter lens A 1835 (z=0.25) spanning 0.35<A O bs<5/im. Sta- 
tistically significant flux is not detected in any of these 
data, including our independent re-analysis of P04's dis- 
covery H- and /T-band data. The 3er detection thresh- 
olds are: F(A o b s =0.44/im)>4.5xl0~ 19 ergs" 1 cm" 2 per spec- 
tral resolution element, R 1W >21 .0, H>25.0, K>25.0, 
F(3.6^m)<0.75/iJy and F(4.5/im)<0.75/iJy, where the pho- 
tometric limits are calculated in apertures with a diameter 3 x 
that of the seeing disk. 

Combining these results with those of Bremer et al. (2004) 
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and Weatherley et al. (2004), we are therefore forced by the 
data to conclude that there is no statistically sound evidence 
for the existence of #1916. We also show that inefficient bad- 
pixel rejection and issues relating to the calculation of the 
background noise can broadly account for the differences be- 
tween P04's near-infrared photometry and our own using the 
same data. The balance of probability is therefore that #1916 
was a false detection in P04's analysis. From a broader per- 
spective, the demise of #1916 warns of the hazards of oper- 
ating close to the detection threshold of deep ground-based 
near-infrared data. 

The need for gravitational magnification to boost the ob- 
served flux of faint (L^O.IL*) galaxies at z^6-8 and popula- 
tions of galaxies at still higher redshifts (see ^ is undimin- 
ished by our results on #1916. However an important issue 
is whether it is feasible to find such galaxies using ground- 
based facilities. We explore this issue using HYPERZ to fit 
synthetic spectral templates to representative data. Our main 
conclusions are that deep near-infrared imaging similar to that 
presented by P04 in combination with a single sensitive op- 
tical non-detection from HST imaging and moderate depth 
Spitzer/IRAC imaging at 3.6 and 4.5/im can discriminate be- 
tween low (e.g. z^2-3) and high (e.g. z>7) redshift solu- 
tions with strong statistical significance. In summary, ground- 
based near-infrared surveys of massive galaxy cluster lenses 



with 10-m class telescopes remain a powerful tool for the dis- 
covery of intrinsically faint galaxies (L<0.1L*) at z>7 that 
may be responsible for cosmic re-ionization. Future surveys 
should combine these data with sensitive space-based optical 
and mid-infrared observations. 
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